function [H0] = Height0(B2bot, k)
%
% Scale Height
%
%DESCRIPTION:
%This function computes the Scale Height.
%
%PROTOTYPE:
% [H0] = Height0(B2bot, k)
%
%--------------------------------------------------------------------------
% INPUTS:
%   B2bot      [1x1]       F2-Layer Bottom Thickness [km]
%   k          [1x1]       Shape Parameter           [-]
%--------------------------------------------------------------------------
% OUTPUTS:
%   H0         [1x1]       Scale Height              [km]
%--------------------------------------------------------------------------
%
%NOTES:
% (none)
%
%CALLED FUNCTIONS:
% (none)
%
%UPDATES:
% (none)
%
%REFERENCES:
% [1] "Ionospheric Correction Algorithm for Galileo Single-Frequency Users"
%      - European GNSS (Galileo) Open Service
% [2] "Electron Density Models and Data for Transionospheric Radio
%      Propagation" - Report ITU-R P.2297-1 (05/2019)
%
%AUTHOR(s):
%Luigi De Maria, Matteo D'Addazio, 2022
%

%% Main Code

%Auxiliary Parameter Ha
Ha = k * B2bot;         %[2] says it is directly H0

%Auxiliary Parameters x and v
x = (Ha - 150) / 100;
v = (0.041163*x - 0.183981)*x + 1.424472;

%Scale Height
H0 = Ha/v;

end